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I. INTRODUCTION 



It is shown that the Collins mechanism explanation of the transverse single- 
Q\ \ spin asymmetries in p^p — > nX leads to a transversely polarized d quark den- 

sity Ard(x) which violates the Soffer bound when one uses several standard 
C**) ' forms for the longitudinally polarized d quark density Ad(x) obtained from 

polarized deep inelastic scattering. Imposition of the Soffer bound with these 
Ad(x) yields results in hopeless disagreement with the data. Remarkably, im- 
position of the Soffer bound, but using parametrizations of Ad(x) that respect 
the PQCD condition Aq(x)/q(x) — * 1 as x — * 1, leads to an excellent fit to 
most of the data. The implications for the polarized DIS neutron longitudinal 
t asymmetry A™ at large x are dramatic. 

PACS numbers: 13.60.Hb, 13.85.Ni, 13.87.Fh, 13.88.+e 

One of the major challenges to the QCD-parton model is the explanation of the large 
(20-40%) single-spin asymmetries found in many semi-inclusive hadron-hadron reactions, of 
which the most dramatic are the polarization of the lambda in pp — > AX and the asymmetry 
under reversal of the transverse spin of the proton in p^p — ► irX. The challenge arises from 
the fact that in the standard approach the basic "two parton — ► two parton" reactions 
involved in the perturbatively treated hard part of the scattering do not possess this kind 
of asymmetry. 

Already some time ago, Efremov and Teryaev jl|] suggested a mechanism for these asym- 
metries utilizing "three parton — > two parton" amplitudes for the hard scattering. This, 
however, necessitates the introduction of a new unknown soft two-parton density, namely 
the correlated probability of finding in the polarized proton a quark with momentum frac- 
tion X\ and a gluon with fraction X2- This quark-gluon correlator contains the dependence 
on the transverse spin of the proton. A fully consistent application of the approach has not 
yet been carried out, though a significant step in the direction has been taken recently by 
Sterman and Qiu ||. 
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Some time ago Sivers || and Collins [Q suggested mechanisms for the asymmetries, which 
are within the framework of the standard "two parton — > two parton" picture, but in which 
the transverse momentum kx of the quark in the hadron plays an essential role. Sivers 
introduces a parton density which depends on the transverse spin St of the proton in the 
form St ■ (kT x p) where p is the momentum of the polarized proton. However such a 
structure is forbidden by time reversal invariance 0. Nonetheless, on the grounds that the 
time-reversal argument could be invalidated by initial or final state interactions, Ansclmino, 
Boglione and Murgia ^,^| have applied the Sivers mechanism to p^p — > irX and shown that 
a very good fit to the data at large xf can be achieved. The problem is that this approach 
raises a fundamental question: if the initial and final state interactions are important then 
it is hard to see why the underlying factorization into hard and soft parts is valid. (The ad- 
vantage of the quark-gluon correlator approach is that it effectively provides initial and final 
state interactions calculable at the parton level). In the Collins mechanism the asymmetry 
arises in the fragmentation of a polarized quark into a pion with transverse momentum pt, 
described by a fragmentation function of the form A N D(z,pr) (see Section II) which is 
convoluted with the transverse spin dependent parton density Axq, about which nothing is 
known experimentally. This structure is not in conflict with time reversal invariance. The 
definition of transverse spin dependent parton densities and possible methods of determining 
them have been studied by Artru and Mekhfi Q, Cortes, Pire and Ralstone |8| and Jaffe 
and Ji Note that in the latter they are referred to as " transversity distributions" and 
that notation differs amongst all these papers. 

An estimate of the size of the Collins effect was first made by Artru, Czyzewski and 
Yabuki |To|], but more recently Anselmino, Boglione and Murgia f| have demonstrated that 
an excellent fit to the data on p^p — » irX can be obtained with the Collins mechanism. 
However their fit is problematic since the A-rg's used in the fit violate the Soffer bound JTT| ] 

|A T g(z)l < ~[q(x) + Aq(x)] (1) 

at large x when Aq(x), the usual longitudinal polarized parton density, is taken from any 
of the standard parametrizations of the longitudinal and unpolarized densities: 

• GRSV-GRV = Gliick, Reya, Stratmann and Vogelsang fjj] + Ghick, Rcya, Vogt @, 

• GS-GRV = Gehrmann and Stirling @ + Gliick, Reya, Vogt ||, 

• LSS-MRST = Leader, Sidorov and Stamenov jl5| + Martin, Roberts, Stirling and 
Thorne 0. 

Note that in the above parametrizations the polarized densities are linked to particular 
parametrizations of the unpolarized densities, as indicated. 

The key point is that the tt~ data demand a large magnitude for Aid at large x, whereas 
Ad(x) is almost universally taken negative for all x, thereby making the Soffer bound much 
more restrictive for the d than for the u quark. 

This raises an intriguing question. There is an old perturbative QCD argument jlTj that 
strictly in the limit x — ► 1 

^-1, (2) 

which would imply that Ad(x) has to change sign and become positive at large x. The 
polarized DIS data certainly require a negative Ad(x) in the range 0.004 < x < 0.75 where 
data exist, but there is no reason why Ad(x) should not change sign near or beyond 0.75. In- 
deed there is a parameterization of the parton densities by Brodsky, Burkhardt and Schmidt 
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[BBS] pl| which has this feature built into it. The original BBS fit is not really competitive 
since evolution in Q 2 was not taken into account, but a proper QCD fit based on the BBS 
parameterization was shown by Leader, Sidorov and Stamenov Jl9| to give an adequate fit 
to the polarized DIS data. 

In this paper we address the question of the correct use of the Collins mechanism in 
which the Soffer bound is respected. We find that it is impossible to get a good fit to 
the tt^~ data when the magnitude of Axd is controlled by (1) in which Ad(x) from any of 
the standard forms given above is used. On the contrary, and most surprisingly, we find 
that parametrizations in which Ad{x)/d{x) — > 1 as x — ► 1 allow a Axd(x) that leads to an 
excellent fit to most of the pion data. 

In Section II we briefly describe the Collins mechanism, and present our results in Section 
III. Conclusions follow in Section IV. 



II. THE MODEL 



As mentioned in the Introduction, we require the Soffer inequality, Eq.(^j), to be re- 
spected by the distribution functions Atu and Axd determined by our fit. Besides, the 
positivity constraint A N D(z) < 2D(z) must hold, since D(z) = \ [D^ (z) + D*{z)] and 
A N D(z) = [D'(z) — D*(z)]. Therefore, the parametrizations are set so that these condi- 
tions are automatically fulfilled, in the following way. First we build a simple function of 
the form x a (l — x) b (or z a {\ — Zf as appropriate), where the powers a, b or a, j3 are > 0, 
and we divide by their maximum value. By allowing an extra multiplicative constant factor 
to vary from —1 to 1, we obtain a function which, in modulus, will never be larger than 1. 
Then we parameterize Axq(x) and A N D(z) by multiplying the functions we built by the 
constraint given by the Soffer inequality or the positivity limit. In this way we make sure 
that the bounds are never broken. 

For the transversity distribution functions we set 

A T u(x) = N u X a{l J bl X) b [\ [u{x] + Au(x)]}, (3) 



A T d(x) = N d 



x c (l 



(c+dy 



) -{\[d{x)+Ad{x)]}.^ 



(4) 



with \N Ui d\ < 1. Here q{x) and Aq{x) are the whole distribution functions, i.e. they contain 
valence and sea contributions (but this is irrelevant at large x since there the contribution 
of the sea is negligible). As in the previous calculation, only u and d contributions are taken 
into account in the polarized proton, so that 



A T u(x) = A T d(x) = A T s(x) = A T s(x) = 0. 



(5) 



For the functions q(x) and Aq(x) we use, for comparison, the "standard" parton 
parametrizations mentioned in Section I and two further parametrizations, one due to Brod- 
sky, Burkhardt and Schmidt (BBS) |l8| which ignores Q 2 -evolution ; and a more consistent 
version of this, due to Leader, Sidorov and Stamenov (LSS)bss (01 which includes the 
Q 2 -evolution. These will be explained in more detail in Section III. 

For the fragmentation function we have 



A N D(z) = N F 



z a (l-z) - 

a a f3' 3 _ 

(Q+/3)° + f j 



2D(z) 



(0) 
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with \Np\ < 1. Here we take into account only valence contributions, so that isospin 
symmetry and charge conjugation give 

A N D%+ = A N D^- = A N D(z) , (7) 



A N D*+ = A N D^- = (8) 

and 

A N D% = A N D$ = A N D% = A N Dio = 1 A w £(z). (9) 

Notice that A w £> is, in fact, a function of the intrinsic transverse momentum p T . 
A N D h / a r (z, p T ) is defined as the difference between the number density of hadrons h, a 
pion in our case, with longitudinal momentum fraction z and transverse momentum p T , 
originating from a transversely polarized parton a with spin either f or j, respectively 

& N D h/a i (z, p T ) = D h/a i {z,p T )~ D h/a i (z,p T ) 

= D h/a i (z,p T )- D h/a i (z, -p T ) , (10) 

where the second line follows from the first one by rotational invariance. Details on the 
integration over the transverse momentum p T and its dependence on z are given in Ref. j(| 
(see Eqs. (17) and (19)). The unpolarized fragmentation function for pions, D(z), is taken 
from Binnewies et al. [M . 

With these ingredients we are now ready to calculate, in complete analogy with Ref. ||, 
the p^p — ► irX single spin asymmetry 

A — ^ ~ ^ (11) 



Here 



J 7TZ 

a,b,c : d 

x A T q a (x a ) q\x b ) A NN a ab ^ cd (x a , x b ,p T ) A N D„ /c (z, p T ) , (12) 



where 



and 



A^ro)b — ^c 1 d J~a^b — >c d 
a -ab^cd /i o\ 

A NN a = — — , (13) 

dt dt 



dal + da i = 2 da un P = 2 ]T / (14) 

x q a (x a ) q b (x b ) -» {x a ,x b ) D w / C (z) ■ 

dt 

All details about the calculation can be found in Ref. J(|. The relation between the above 
notation and that of Q is : q a {x) = f a / P (x) and At9°(x) = P a / pT f a j p (x). 
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III. RESULTS 



We start by running two fits to the E70A experimental data |2lJ using the popular GS 
Jl4[ polarized densities in conjunction with the GRV (T^| unpolarized densities, and the 
latest LSS |l5|] polarized densities in conjunction with the MRST jl6| unpolarized densities. 
It should be noted that the 1996 analysis of Gehrmann and Stirling was done prior to 
the publication of a great deal of new, high precision data on polarized DIS, whereas the 
Leader, Sidorov and Stamenov analysis includes all of the present world data. Fig. 1 shows 
the complete failure to fit the data when the Soffer bound is implemented using the GS-GRV 
parametrizations (\ 2 /DOF = 25 !). The corresponding transverse densities Atu(x) and 
AtcI(x) are shown in Fig. 2. In this fit only \Np\ = 1 is possible and Fig. 2 corresponds to 
Np — — 1. The sign is discussed later. 

A somewhat better picture emerges when using the LSS-MRST results to implement the 
Soffer bound: Fig. 3. The fit looks reasonable out to sj? = 0.4 but fails beyond that. One 
finds x 2 /DOF = 6.12. The transverse densities are shown in Fig. 4, where the curves 
correspond to negative Np. 

Note that in both these fits one finds a = (3 = in Eqn. (6), showing that the magnitude 
of A N D{z) is maximized at each z-value. 

The reason for the failure of GS-GRV case and for the relative success of LSS-MRST 
can be understood by observing in Fig. 5 that the Soffer bound on Ayd(x) is much more 
restrictive at large x in the GS-GRV case. Comparison of Fig. 5 with Fig. 6 also indicates 
the source of the problem. The asymmetries for ir^ are of roughly equal magnitude whereas 
the Soffer bound restrictions are much more severe for the d quark as a consequence of Ad(x) 
being negative for all x. 

This suggests an intriguing possibility. The polarized DIS data only exist for x < 0.75 
and there is really very little constraint from DIS on the Aq(x) for x near to and beyond 
this value. At the same time there are perturbative QCD arguments [n7| which suggest that 



Aq(x) 



1 as x — > 1 (15) 



q{x) 

and, indeed, even more precisely, that 

q(x) - Aq(x) cx (1 — x) 2 q(x) as x — *• 1 . (16) 

This constraint is almost universally ignored in parameterizing the Aq(x), on the grounds 
that (|l^) is incompatible with the evolution equations. But this is a "red herring" since the 
evolution equations do not hold in the region where (|l6| ) is valid, approaching the border of 
the exclusive region. 

The imposition of ( |l5| ) is exactly what we need for Aq(x) to change sign and become 
positive at large x, thereby diminishing the restrictive power of the Soffer bound on Axdjx). 
In fact there does exist a parametrizations of the Aq(x) which respects (|l5| ) and (16) 



namely that of Brodsky, Burkhardt and Schmidt (BBS) [T8|. Unfortunately BBS did not 
include any (Revolution where determining the numerical values of their parameters from 
the DIS data, so their fit is not really adequate. 

However, Leader, Sidorov and Stamenov |l9[ ] made an extensive study, using the BBS 
functional forms, but including Q 2 -evolution, and found a very good fit (LSS)bbs to the 
polarized DIS data available in 1997. 

It can be seen in Fig. 5 that the Soffer bound on Ard(x) is much less restrictive for the 
BBS case and that the (LSS)fjgg bound is rather similar to that of the LSS-MRST case, but 
is less restrictive for x > 0.7. It is important to realize that although the Axd(x a ) needed 
in ( |l2] ) are tiny for such large values of x a , this is compensated for by the fact that large Xf 
then demands very small Xb, where the unpolarized densities grow very large. 
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GS-GRV 



LSS-MRST 



BBS 



(LSS)bbs 



N F N U 
N F N d 



-0.43 
1.00 
4.33 
0.00 
0.00 
0.00 
25 



-0.73 
1.00 
3.03 
0.00 
3.48 
0.76 
6.12 



0.54 
0.88 
3.17 
0.00 
3.57 
0.00 
1.45 



-0.49 
0.91 
3.46 
0.00 
3.32 
0.00 
2.41 



a 



b 



c 



d 

X 2 /DOF 



TABLE I. Parameters determined by the fit in the four different parameterization schemes and 
the corresponding x 2 /DOF. 

Note that the relative signs of N u and Nd are opposite, but their absolute signs are not 
determined since, in principle, Np can be positive or negative. However, if one uses an 
SU(6)f wave function for the proton, one finds Atu positive and A^e? negative, so it seems 
reasonable to hypothesize that N u > and Nd < 0. For this reason we have chosen Np to 
be negative in the above. Note that N u and Nd are not a direct measure of the magnitudes 
of Atu and Axd. Their role is linked specifically to the Soffer bound. The relative behavior 
of Atu and Ard can be seen in Figs. 2, 4, 9, 10. 

Indeed, as expected, we find that a significantly better fit to the asymmetry data is 
achieved using the BBS and the (LSS)bbs parametrizations, with \ 2 /DOF = 1.45 and 
X 2 /DOF = 2.41 respectively. As can be seen in Fig. 7 and 8. The curves reproduce the 
trends in the data right out to sf ~ 0.7. Figs. 9 and 10 show how similar the allowed ranges 
of transverse polarized densities are in the two cases. In Fig. 9, 0.88 < \Np\ < 1 , whereas 
in Fig. 10 0.91 < \Np\ < 1 . As before the curves correspond to negative Np. 



The parameter values for all the parametrizations are shown in Table 1, where it should 
be recalled \N F \, \N U \ and |AT d | < 1. 



We have demonstrated that the Collins mechanism is able to explain much of the data 
on the transverse single spin asymmetries in p^p — > ttX, namely the data in the region 
xp < 0.7, if, and only if, the longitudinal polarized c?-quark density, which is negative for 
small and moderate x, changes sign and becomes positive at large x. There is hopeless 
disagreement when using the longitudinal polarized densities due to Gehrmann and Stirling 
pi , and matters are significantly better when using the most up to date parameterization of 
Leader, Sidorov and Stamenov |15| . But the most successful fits arise from parametrizations 
Ji]|[l9| which respect the PQCD condition Ad(x)/d(x) ->■ 1 as x — ► 1. 

For parametrizations of Ad(x) with this property there are interesting consequences in 
polarized DIS, namely, the neutron longitudinal asymmetry A™(x) should change sign and 
tend to 1 as x — > 1 (see Fig. 11). The region of large x has hardly been explored in polarized 
DIS up to the present. Clearly a study of this region might turn out to be very interesting. 

There remains, however, the problem of the p^p — > irX data at the largest values of 
Xp so far measured, i.e. 0.7 < Xp < 0.82. It does not seem possible to account for 
these asymmetries within the framework of the Collins mechanism. On the other hand Qiu 
and Sterman pi, using a "three parton — » two parton" amplitude for the hard partonic 
scattering and a "gluonic pole" mechanism, claim to be able to reproduce the very large 
asymmetries at ij? ~ 0.8. However their study must be considered as preliminary, since 



IV. CONCLUSIONS 
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it relies on a completely ad hoc assumption that the essential new twist-three quark-gluon- 
quark correlator function Tp\x,x), for given flavour /, is proportional to qf(x), and no 
attempt is made to fit the detailed x ^-dependence of the data. 
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FIG. 1. Single spin asymmetry for pion production in the process p^p — > ttX as a function of 
xf, obtained by using the GS-GRV |l3||l4[ sets of distribution functions. The solid line refers to 
7r + , the dashed line to 7r° and the dash-dotted line to 7r~. 
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FIG. 2. The distribution functions Atu(x) and At<1(x), as obtained by using the GS-GRV 
[ p^|Jl4| ] distribution functions. The dotted lines are the boundaries imposed by the Soffer inequality. 
For Ard(x) the dotted line is invisible since Ard(x) completely saturates the Soffer bound in the 
whole x region. For discussion of signs, see text. 
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FIG. 3. The single spin asymmetry for pion production in the process p^p — > tyX as a function 
of xp, obtained using the LSS-MRST Jl5] , [l6| distribution functions. The solid line refers to ir + , the 
dashed line to n° and the dash-dotted line to tv~ . 
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FIG. 4. The distribution functions Atu(x) and Ard(x) versus x, as determined by the fit using 
the LSS-MRST ( fl5||l6| ] distribution functions. The dotted lines are the boundaries imposed by the 
Soffer inequality. For discussion of signs see text. 
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FIG. 5. The boundaries imposed on Ard by the Soffer inequality. Note that the GS-GRV 
distribution functions give a much tighter bound than any of the other parameterization sets. 
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FIG. 6. The boundaries imposed on Atii by the Soffer inequality. For the u flavour the Soffer 
bound is very similar in each set of parametrizations. 



10 




-0.6 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

X F 

FIG. 7. The single spin asymmetry for pion production in the process p^p — * tyX as a function 
oi xf, obtained by using the BBS |l8) distribution functions. The solid line refers to n + , the dashed 
line to 7T° and the dash-dotted line to tt~ . 
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FIG. 8. The single spin asymmetry for pion production in the process p^p — > tyX as a function 
of xf, determined by the fit using the (LSS)sss Jit]] distribution functions. The solid line refers to 
7r + , the dashed line to ir° and the dash-dotted line to ir~ . 
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FIG. 10. The range of allowed distribution functions Atu(x) and Ard(x) versus x, as deter- 
mined by the fit using the (LSS)bbs (is| distribution functions. The dotted lines are the boundaries 
imposed by the Soffer inequality. For signs see discussion in text. 
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FIG. 11. The neutron longitudinal asymmetry A"(x) for Q 2 ~ 1 — 4 GeV 2 , as obtained by using 
the BBS and (LSS)sss parametrizations (solid and dashed lines respectively), and the LSS-MRST 
parametrizations (dash-dotted line). 
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